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Abstract — The  formation  of  hafnium  oxide  (Hf02)  layers  in  slit  structures  by  pulsed  metalorganic  chemical- 
vapor  deposition  (MOCVD)  method  with  discrete  dosage  of  reactants  has  been  studied.  The  MOCVD  pro¬ 
cess  employed  hafnium  tetrakisdipivaloylmethanate  as  a  precursor  and  oxygen  as  a  gaseous  reactant.  The  slit 
structure  was  formed  by  a  silicon  plate  and  a  glass  substrate  with  a  chromium- film  pattern.  The  slit  width  was 
varied  within  0. 1—2.0  mm,  so  that  the  aspect  ratio  changed  from  150  to  30.  The  deposited  layers  were  studied 
by  the  scanning  ellipsometry  technique.  The  phenomenon  of  image  transfer  from  the  glass  substrate  to  the 
opposite  plate  in  the  form  of  a  nanosized  Hf02  film  relief,  which  reproduces  the  chromium-film  pattern,  is 
discovered.  The  in-plane  resolution  of  the  obtained  image  in  these  experiments  reached  0.025  mm.  This 
effect  can  be  used  for  nonlithographic  formation  of  images  of  film  structures. 
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Artificial  ordered  microsystems — e.g.,  of  the  types 
created  by  MEMS  and  LIGA  technologies,  etc. — are 
finding  increasing  use  in  modem  advanced  industries. 
In  some  cases,  the  fabrication  of  these  microsystems 
requires  applying  functional  nanodimensional  coat¬ 
ings  onto  a  developed  microrelief  containing  slits, 
channels,  and  wells  with  high  aspect  (width/length) 
ratios  up  to  1  :  1000  [1,  2], 

The  present  work  was  aimed  at  studying  specific 
features  of  the  deposition  of  nanosized  metal  oxide 
layers  in  a  slit  system  with  a  high  aspect  ratio. 

Layers  of  hafnium  oxide  (Hf02)  were  formed  by 
the  method  of  pulsed  metalorganic  chemical -vapor 
deposition  (MOCVD)  employing  discrete  dosage  of 
reactants  without  a  carrier  gas.  Volatile  hafnium  tet¬ 
rakisdipivaloylmethanate  Hf(thd)4  (where  thd  = 
(CH3)3CCOCHCOCC(CH3)3)  was  used  as  a  meta¬ 
lorganic  precursor  and  oxygen  was  a  gaseous  reactant. 
The  MOCVD  apparatus  has  been  described  in  detail 
elsewhere  [3].  The  substrate  temperature  during 
hafnium-oxide-layer  growth  was  maintained  at  750  ± 
2  K,  while  the  Hf(thd)4  vapor-source  temperature  was 
500  ±  2  K,  which  corresponded  to  a  partial  vapor  pres¬ 
sure  of  212.8  Pa.  The  oxide-layer  growth  was  carried 
out  under  conditions  of  discrete  oxygen  dosage  by 
admission  to  an  initial  pressure  of  2  x  103  Pa  in  the 
reactor  chamber.  The  residual  gas  pressure  in  the  reac¬ 
tor  upon  every  cycle  was  1.33  Pa.  The  total  number  of 


deposition  cycles  was  the  same  (100)  in  all  experi¬ 
ments. 

The  experiments  with  Hf02  film  deposition  were 
performed  using  a  slit  structure  formed  by  two  parallel 
plates  with  a  controlled  gap  width  (Fig.  la).  The  upper 
plate  comprised  a  2-mm-thick  glass  substrate  with 
dimensions  of  30  x  30  mm  bearing  a  pattern  of  chro¬ 
mium  (Cr)  thin-film  squares  (Figs,  la,  lb).  The  Cr 
film  thickness  was  25  nm.  The  lower  plate  was  made  of 
a  Si(100)  wafer.  The  slit  width  in  experiments  could  be 
varied  within  0.2— 1.0  mm,  so  that  the  aspect  ratio 
changed  from  150  to  30.  During  the  deposition  of 
oxide  layers,  the  Si  substrate  could  be  heated  via  con¬ 
tact  with  a  heater. 

The  thickness  of  hafnium  oxide  films  deposited 
onto  silicon  substrates  was  determined  by  the  method 
of  laser  ellipsometry  on  the  MIKROSKAN  high-reso¬ 
lution  scanning  ellipsometer  developed  at  the  Institute 
of  Semiconductor  Physics  (Novosibirsk).  The  ellipso- 
metric  angles  were  measured  at  a  light  wavelength  of 
632.8  nm  for  a  beam-incidence  angle  of  60°.  The  spa¬ 
tial  resolution  in  the  scanning  ellipsometry  regime  was 
about  10  pm. 

During  the  formation  of  a  hafnium-oxide  layer  on 
the  lower  (Si)  substrate,  the  deposit  exhibited  a  visible 
image  that  corresponded  to  the  pattern  of  chromium 
film  on  the  upper  (glass)  substrate  (Fig.  lc).  The  thick¬ 
ness  of  the  Hf02  layer  formed  against  the  Cr  squares 
was  smaller  than  that  against  the  gaps  between  squares 
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Fig.  1.  Patterned  oxide-layer  deposition  in  the  MOCVD  system:  (a)  slit  structure  formed  by  (7)  glass  plate  and  (2)  silicon  sub¬ 
strate,  (b)  fragment  of  glass  plate  with  lithographic  pattern  of  Cr  thin-film  squares,  and  (c)  photograph  of  the  image  transferred 
to  hafnium-oxide-covered  Si  substrate. 


(Fig.  2).  We  have  also  studied  the  influence  of  the  slit 
width  (0.2,  0.5,  and  1  mm)  on  the  contrast  of  trans¬ 
ferred  images.  As  the  slit  width  was  increased  to  1  mm, 
the  relief  of  the  deposited  film  exhibited  smearing.  The 
oxide  layer  was  also  formed  on  the  upper  (glass)  plate 
with  Cr  squares,  but  this  deposit  could  be  studied  by 
laser  interferometry  in  regions  free  of  the  metal 
because  the  refractive  indices  of  the  glass  substrate  and 
deposit  were  close.  At  a  slit  width  of  0.5  mm,  the 
obtained  image  resolution  reached  0.1  mm. 

Thus,  our  experiments  on  pulsed  MOCVD  with 
discrete  dosage  of  reactants  in  a  slit  structure  revealed 
the  phenomenon  of  image  transfer  from  one  patterned 
plate  to  a  parallel  substrate. 

In  order  to  explain  the  observed  phenomenon,  we 
took  into  account  the  following  peculiarities  of  the 
experiment. 

(1)  The  observed  image  is  formed  as  a  result  of  the 
surface  chemical  reaction  of  decomposition  of  a  met- 
alorganic  compound.  This  surface  reaction  (with 
allowance  for  heat  supply  from  the  substrate  to  the  gas¬ 
eous  phase)  is  accompanied  by  the  active  generation  of 
new  components  in  the  vapor— gas  mixture,  which 
leads  to  pressure  variations  in  the  near-surface  layer 
and  the  related  flux  of  reaction  products  from  the  sur¬ 
face  to  space  between  plates. 


(2)  The  values  of  thermodynamic  parameters  in  the 
system  are  such  that  the  mean  free  path  of  molecules 
is  about  two  orders  of  magnitude  smaller  than  the  slit 
width  (i.e.,  the  Knudsen  number  is  on  the  level  of 
0.01). 

(3)  The  relationship  between  the  pressure  in  the 
reactor  and  the  slit  size  is  such  that  the  time  of  gas  flow 
into  the  slit  is  much  shorter  than  the  period  of  expo¬ 
sure  during  one  MOCVD  cycle. 

The  above  limitations  do  not  allow  the  observed 
synergetic  effect  of  image  transfer  to  be  considered 
using  standard  approaches  and  approximations  such 
as  the  models  of  (i)  ultrararefied  gas,  (ii)  continuous 
medium,  and  (iii)  convective  flows  and  the  theory  of 
Benard  convection  instability.  In  view  of  these  limita¬ 
tions,  it  was  necessary  to  carry  out  theoretical  analysis 
of  the  observed  physical  phenomenon  using  a  mathe¬ 
matical  model  with  a  minimum  set  of  phenomenolog¬ 
ical  assumptions.  It  was  found  that  the  most  accept¬ 
able  approach  was  offered  by  the  method  of  molecular 
dynamics  with  the  potentials  of  intermolecular  inter¬ 
action  and  molecule— substrate  surface  interaction 
determined  by  quantum-mechanical  calculations. 

The  system  geometry  was  modeled  by  two  coaxial 
disks  of  the  same  radius  R2,  occurring  in  a  cylindrical 
reactor  of  radius  R{  (R2  <  R\)  and  height  H.  The  inter- 
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Fig.  2.  Relief  of  Hf02  film  on  Si  substrate  by  data  of  scan¬ 
ning  ellipsometiy:  (a)  flat  image  (with  gray-tint  scale  of 
heights);  (b)  three-dimensional  relief. 


action  of  molecules  was  described  using  a  model  of 
spheres  wit  radius  r,  interacting  according  to  the  Len- 
nard- Jones  potential  with  parameters  du  =  3.405  A 
and  eLJ  =  1.65  x  10-14  erg.  The  interaction  of  mole¬ 
cules  with  the  surfaces  of  disks  and  reactor  walls  was 
described  by  the  Lennard- Jones  potential  branch  cor¬ 
responding  to  repulsion. 

The  flow  of  molecules  in  the  slit  at  constant  gasdy- 
namic  parameters  (corresponding  to  experimental 
conditions)  in  the  reactor  was  simulated  as  follows.  At 
the  initial  moment  of  time,  the  positions  and  velocities 
(momenta)  of  molecules  in  the  ring-shaped  space 
from  R  |  to  R2  were  set  with  the  aid  of  a  random  number 
generator.  The  spatial  arrangement  of  molecules  was 
uniform,  while  the  momenta  were  distributed  accord¬ 
ing  to  the  Maxwell  law  for  gas  temperature  Tg.  The 
number  of  molecules  was  determined  in  accordance 
with  a  simple  ideal  gas  model.  The  system  was  allowed 
to  evolve  by  itself  with  the  molecules  moving  accord¬ 
ing  to  the  laws  of  classical  mechanics.  Molecular  tra¬ 
jectories  in  time  were  calculated  using  the  classical 
equations  of  motion.  As  soon  as  a  molecule  fell  in  the 
space  between  disks  (i.e.,  was  lost  from  the  reactor), 
new  molecules  were  injected  into  the  reactor  with  the 
aid  of  a  random  number  generator  so  as  to  satisfy  two 
conditions:  (i)  the  total  number  of  molecules  in  the 
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Fig.  3.  Temporal  variation  of  heat  flux  Q  to  the  center  of  the 
lower  disk  calculated  for  two  sets  of  the  probabilities  of  heat 
transfer  to  the  surface:  (i)  JVe  =  0. 1 ,  Wt  =  0.9  and  (2)  We  = 
0.99,  Wt  =  0.01. 


system  must  be  below  theoretical  and  (ii)  the  number 
of  molecules  in  the  reactor  must  be  below  theoretical. 
If  the  total  number  of  molecules  in  the  system  is  equal 
to  theoretical,  the  injection  ceases  and  the  number  of 
molecules  in  the  reactor  and  between  disks  exhibits 
fluctuational  variation.  The  energy  exchange  with  the 
surface  of  the  reactor  and  disks  was  simulated  using  the 
method  described  in  [4] .  The  main  parameters  of  this 
energy  exchange  include  the  probability  of  elastic 
scattering  on  the  surface  ( We)  and  the  probability  of 
scattering  with  heat  transfer  (Wt).  In  the  absence  of 
chemical  reactions,  these  quantities  must  obey  the 
condition  We+  Wt=  1, 

Analysis  of  the  gas  concentration  and  temperature 
between  disks  was  performed  as  follows.  A  set  of  disk 
radii  was  introduced  such  that  the  concentric  rings 
between  these  radii  had  equal  areas  (and  the  corre¬ 
sponding  mesocells  between  disks  had  equal  volumes). 
Then,  the  state  of  gas  in  these  mesocells  was  analyzed 
at  certain  time  steps  by  determining  the  numbers  of 
molecules  in  these  mesoscopic  volumes,  their  temper¬ 
atures,  and  the  fluxes  of  molecules  and  energy  to  the 
surface  of  mesoscopic  rings  on  the  upper  and  lower 
disks.  In  addition,  the  fluxes  of  momentum  compo¬ 
nent  perpendicular  to  the  surface  was  determined, 
which  allowed  the  pressure  on  the  surface  of  the  reac¬ 
tor  and  mesoscopic  rings  to  be  calculated. 

Figure  3  shows  typical  results  of  calculations  of  the 
heat  flux  to  mesoscopic  rings  on  the  upper  and  lower 
disks  for  two  sets  of  the  probabilities  of  heat  transfer  to 
the  surface:  (1)  We  =  0. 1  and  Wt  =  0.9  (for  the  entire 
surface  of  the  reactor  and  disks)  and  (2)  We  =  0.01  and 
Wt  =  0.99  (in  the  central  part  of  the  upper  disk,  which 
models  the  properties  of  Cr  coating)  and  We  =  0.99 
and  Wx  =  0.01  on  the  rest  of  the  surface  (and  the  prob- 
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ability  of  heat  exchange  in  the  reactor  being  the  same 
as  in  the  first  case).  The  distribution  of  temperatures  in 
the  system  was  as  described  above. 

The  results  of  model  calculations  showed  that  the 
probabilities  of  heat  exchange  between  molecules  and 
the  surface  influence  the  heat  flux  to  mesoscopic  rings 
on  the  lower  disk,  which,  in  turn,  affects  the  rate  of 
surface  reactions  and,  hence,  the  amount  of  a  depos¬ 
ited  reaction  product,  thus  determining  the  image 
transfer  from  the  upper  to  lower  disk. 

Thus,  it  was  established  that  the  intense  heat  and 
mass  transfer  under  strongly  nonequilibrium  condi¬ 
tions  in  the  gas  phase  during  a  pulsed  MOCVD  pro¬ 
cess  in  slit  structures  results  in  the  phenomenon  of  a 
synergetic  transfer  of  images  from  one  plate  to  the 
opposite  substrate,  which  proceeds  according  to  a  pre¬ 
viously  unknown  mechanism.  Direct  modeling  by  the 
method  of  molecular  dynamics  revealed  an  asymme¬ 
try  of  heat  fluxes  in  the  system  that  depends  on  the 
thermal  conductivity  and  temperature  diffusivity  of 
the  upper  disk  material,  which  is  evidence  of  the  valid¬ 
ity  of  the  adopted  theoretical  approach  to  explanation 


of  the  observed  phenomenon.  The  obtained  results 
determine  the  direction  of  subsequent  simulations, 
which  must  take  into  account  chemical  reactions 
involved  in  oxide-layer  formation. 
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